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INTRODUCTION 


The  effects  of  weather  on  the  planning  for  and  execution  of  combat  missions 
have  been  apparent  to  military  men  as  long  as  man  has  fought  man.  Histor¬ 
ically,  the  weather  has  been  considered  to  a  great  extent  to  be  a  natural 
obstacle,  akin  to  difficult  terrain,  about  which  soldiers  could  do  little 
but  accept  the  consequences.  And  unlike  natural  obstacles,  the  weather 
often  varied  considerably  and  was  difficult  to  predict;  one  had  to  expect 
and  prepare  for  the  worst. 

Trafficabil ity  is  of  critical  concern  to  the  army  in  the  field.  Battle¬ 
field  commanders  must  have  information  on  fluctuations  in  soil  moisture 
content  to  make  decisions  concerning  tactical  maneuvers  on  the  battlefield. 
Soil  moisture  conditions  at  a  given  time  and  place  are  a  function  of: 

(1)  water  input  to  the  ground  surface  from  precipitation,  (2)  water  loss 
from  the  surface  by  evaporation  and  transpiration,  and  (3)  changes  in 
soil  moisture  storage. 

In  this  study  the  primary  concern  is  with  item  (1)  of  the  preceding  para¬ 
graph.  Knowledge  of  where  and  when  rain  is  occurring  as  well  as  the  rate  of 
rainfall  will  provide  critical  data  on  changes  in  soil  moisture  and  thus  on 
trafficabil ity. 

Much  of  the  research  work  which  is  aimed  at  using  satellite  data  to  arrive 
at  estimates  of  rainfall  amounts  has  been  very  subjective  in  nature;  that 
is,  each  person  viewing  the  satellite  imagery  may  interpret  it  in  a  slightly 
different  manner,  thus  arriving  at  variations  in  the  estimated  surface 
rainfall  amounts.  The  elimination  of  as  much  of  this  subjectivity  as 
possible  is  desirable.  An  example  of  this  subjective  approach  is  Scofield 
and  Oliver's1  work  which  made  use  of  infrared  ( I R )  data  to  make  quanti¬ 
tative  estimates  of  cloud  heights.  Rainfall  estimates  were  made  from 
enhanced  IR  satellite  images  based  on  the  premise  that  the  higher  the 
cumulonimbus  cloud  top  the  heavier  the  rain.  Additional  information  on 
the  rate  of  anvil  growth,  the  position  of  the  cumulonimbus  under  the 
spreading  anvil,  merging  cells,  merging  convective  cloud  lines,  and  over¬ 
shooting  tops  has  been  used  to  improve  the  model.  These  factors  are  then 
incorporated  into  a  decision  tree  to  aid  the  user  in  determining  where 
and  how  much  rain  has  fallen. 

This  study  will  attempt  to  show  that,  by  fitting  particular  nonlinear  auto¬ 
covariance  and  cross-covariance  functions  to  space-time  covariance  values 
calculated  from  satellite  data,  storm  characteristics  such  as  size,  ellip- 
ticity,  motion,  growth,  and  decay  of  a  storm  can  be  determined.  These 
modeled  characteristics  then  can  be  related  to  storm  features  found  by 
Scofield  and  Oliver.1  The  storm  features  in  turn  are  related  to  actual 
rainfall  amounts.  The  results  of  this  study  are  a  first  step  in  removing 
the  subjectivity  from  current  rairfall  estimation  techniques. 


JR.  A.  Scofield  and  V.  J.  Oliver,  1977,  "A  Scheme  for  Estimating  Convec¬ 
tive  Rainfall  from  Satellite  Imagery,"  NOAA  Technical  Memorandum  NESS  86 
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DATA  PROCESSING  PROCEDURES 


The  Geostationary  Operational  Environmental  Satellite  (GOES)  surveys  the 
earth  from  a  nominal  altitude  of  35,800  km.  The  GOES  has  a  rotation  speed 
of  100  rpm  which  generates  a  west-to-east  scan  of  the  earth.  A  stepping 
mirror  allows  the  sensors  to  accomplish  the  north- to-south  scan  of  the 
earth.  Eight  visible  and  two  IR  sensors,  with  resolutions  of  1  by  1  km 
and  4  by  8  km,  respectively,  provide  the  data  base.  In  one  revolution 
of  the  satellite,  the  earth  is  scanned  for  about  18  degrees  with  each  scan 
consisting  of  data  from  one  IR  sensor  and  eight  visual  sensors.  These 
data  are  transmitted  to  the  command  station  at  a  28-megabit/s  rate  and 
then  are  time-stretched  to  fill  the  remaining  342  degrees  of  satellite 
rotation  and  retransmitted  back  to  the  satellite.  GOES  then  retransmits 
the  data  back  to  earth  where  they  are  available  to  the  users  by  a 
2-megabit/s  data  rate. 

When  the  satellite  signal  is  received  at  the  Atmospheric  Sciences  Laboratory 
(ASL)  Direct  Readout  Ground  Station  (DRGS),  it  is  processed  through  the  Bit 
Frame  Synchronizer/Sectionizer  System  (BFSSS).  First  the  signal  is  detected 
and  demodulated  to  digital  baseband  nonreturn  to  zero-level  (NRZL)  code  by 
a  demodulator.  The  signal  is  then  routed  to  two  bit  synchronizers,  one 
dedicated  to  infrared  and  the  other  to  visible.  The  outputs  of  the  bit 
synchronizers,  each  consisting  of  NRZL  data  and  bit  clock,  are  routed  to 
the  digital  interface  where  the  interface  establishes  synchronization  on 
the  IR  signal  and  determines  the  transmission  mode.  Synchronization  is 
then  established  on  the  visible  signal.  The  digital  interface  is  operated 
by  a  PDP  11/45  minicomputer.  Through  a  computer  program,  sectors  of  data 
of  interest  to  a  researcher  are  stored  on  magnetic  tape  or  displayed  in 
real  time. 

Once  the  sector  of  interest  is  determined,  the  IR  and  visible  data  are 
separated  into  two  files  for  further  processing.  Often  when  GOES  data 
are  viewed  on  a  display  monitor,  it  is  obvious  that  entire  lines  of  data 
are  unacceptable  for  use  in  an  analysis  program,  and  random  errors  are 
sometimes  present  in  the  data.  These  errors  require  the  use  of  a  filter¬ 
ing  and  smoothing  program  to  render  the  data  acceptable  for  use  in  analysis 
programs.  In  outline  form  the  filtering  and  smoothing  program  handles 
the  data  in  the  following  manner.  Image  positions  are  defined  by  rectan¬ 
gular  coordinates  (i,  j)  where  i  =  1,  2,  ...,  N  gives  the  relative  hori¬ 
zontal  position  (pixel)  and  j  =  1,  2,  ....  M  gives  the  relative  position 
from  the  top  of  the  image  (line);  and  image  intensity  at  point  (i,  j)  is 
denoted  by  F(i,  j).  Examination  of  the  data  showed  that  bad  isolated 
points  were  characterized  by  F(i,  j)  >  elf  or  F(i,  j)  <  e2.  For  visible 
data,  0!  was  assigned  the  value  59  and  e2  the  value  4;  while  for  IR  data, 

0i  and  02  were  assigned  values  of  252  and  4,  respectively.  An  isolated 
bad  point,  F(i,  j),  is  replaced  by  the  average  of  its  neighbors 
F(i  -  1,  j)  and  F(i  +  1,  j)  on  the  same  line  j. 

Unacceptable  lines  were  characterized  by  rapid  oscillations  of  intensities 
along  a  line  and  are  detected  by  a  large  value  of  the  "average  slope"  of 
the  line  by  equation  (1). 
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F(1  +  1,  j)  -  F(i,  j)  . 


(1) 


The  lines  were  defined  as  bad  when  A  was  greater  than  five.  After  examin¬ 
ing  a  large  number  of  cases,  the  author  found  that  good  lines  have  a  value 
of  A  less  than  three  and  bad  lines  have  a  value  of  A  near  20.  In  actual 
computer  implementation,  lines  are  processed  first  and  random  points  next. 
The  program  also  allows  up  to  three  consecutive  bad  lines  or  bad  points 
before  default  of  the  entire  data  set.  For  the  purposes  of  this  study, 
the  data  were  segmented  so  that  only  one  convective  complex  will  appear 
in  a  data  set.  Tests2  have  shown  that  the  complete  data  set  could  possibly 
show  some  decimation  while  still  producing  the  same  covariance  structure. 
The  results  of  this  study  were  derived  by  using  every  fifth  pixel  and 
every  second  line  in  the  infrared  and  every  twentieth  pixel  and  every 
sixteenth  line  in  the  visible. 

After  this  facet  of  the  data  processing  procedure  was  complete,  latitude 
and  longitude  for  each  point  of  the  decimated  data  set  were  calculated. 

The  data  file  and  the  latitude-longitude  file  were  then  used  to  calculate 
the  raw  covariance  matrices.  A  data  processing  flow  chart  is  shown  in 
figure  1. 


NAVIGATION  FOR  GEOSYNCHRONOUS  METEOROLOGICAL  SATELLITES 

The  navigation  algorithm  used  by  the  ASL  Meteorological  Analysis  System 
(AMAS)  is  a  modification  of  Phillips  and  Smith's  work.3  Similar  software 
packages  are  being  used  by  NASA/GSFC  and  the  Universities  of  Wisconsin  and 
Colorado  State. 

The  purpose  of  the  navigation  software  is  to  provide  the  capability  for 
accurately  calculating  the  latitude  and  longitude  of  any  point  in  the 
image  data  and,  conversely,  for  computing  the  position  on  the  image  of 
any  specified  earth  location.  The  data  required  to  perform  these  trans¬ 
formations  consist  of  four  principal  components:  (1)  spacecraft  orbital 
elements  for  computing  satellite  position  at  any  specified  time,  (2)  space¬ 
craft  spin  axis  attitude,  specified  by  the  right  ascension  and  declination 


2A.  Eddy  and  L.  Hembree,  1977,  "Space-Time  Sampling  from  SMS  Satellite 
Data  Required  to  Define  Convective  Storms,"  Final  Report,  Contract 
DAAG29-76-D-0100,  US  Army  Research  Office 

3D.  R.  Phillips  and  E.  Smith,  1974,  "Geosynchronous  Satellite  Navigation 
Model,"  Internal  Report,  University  of  Wisconsin  Space  Sciences  and 
Engineering  Center 


7 


Figure  1.  Data  processing  flowchart. 
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of  the  spin  vector  in  celestial  coordinates,  (3)  small  angles  of  pitch, 
roll,  and  yaw  that  specify  the  misalignment  between  instrument  optical 
axes  and  spacecraft  spin  axis  and,  (4)  scan  phase  correction  parameters 
for  each  image  to  which  the  navigation  applies.  These  parameters  (item 
(4)  above)  are  the  initial  offset  and  drift  rate,  in  the  pixel  (horizontal) 
direction,  between  the  center  of  the  earth  and  the  center  of  the  master 
image. 

Three  coordinate  systems  are  used  in  the  navigation  solution,  namely: 

(1)  a  celestial  inertial  system  for  description  of  orbit  and  attitude, 

(2)  a  geographic  system  rotating  with  the  earth  for  description  of  land¬ 
marks,  and  (3)  a  time  dependent  satellite  coordinate  system  for  descrip¬ 
tion  of  the  scanning  geometry  (figure  2).  For  a  detailed  description  of 
this  solution,  consult  Billingsley  et  al.4 


THE  OBJECTIVE  ANALYSIS  TECHNIQUE 

The  objective  analysis  technique  used  in  this  project  requires  that  the 
intrarelationships  existing  within  each  single  parameter  data  set,  as 
well  as  the  interrelationships  existing  between  these  data  sets,  be 
determined.  These  relationships  are  then  described  in  terms  of  correla¬ 
tion  and  cross-correlation  functions  which  model  the  spatial  and  temporal 
structure  of  the  phenomenon  as  reflected  in  the  observations  themselves. 

The  analysis  model  used  is  an  extension5  of  classical  multiple  linear 
regression.  This  extension  consists  of  modeling  the  structure  that  is 
representative  of  the  information  contained  in  the  time  series  data  set 
being  analyzed  and  then  using  that  model  in  the  determination  of  the 
regression  weights. 

The  model  for  the  univariate  case  is  given  by 


Y  =  Xe  +  e 


(2) 


where  Y,  the  predictand,  is  an  NX  1  matrix;  X,  the  predictor  set,  is  an 
NXM  matrix;  e,  the  vector  of  regression  weights,  forms  an  MX  1  matrix; 
and  e,  the  noise  vector,  forms  an  NX1  matrix. 

In  classical  regression.  If  E  [e]  =  0,  where  E  is  the  expectation  operator, 
Vcr2  =  E  [cct]  is  the  variance-covariance  matrix  of  e,  where  et  denotes 


4J.  J.  Billingsley  et  al . ,  1977,  "A0IPS  METPAK,  A  Meteorological  Data 
Processing  System,"  Computer  Science  Corporation  Report  CSC/SD-77/6084 

r'A.  Eddy,  1973,  "The  Objective  Analysis  of  Atmospheric  Structure," 

J  Meteorol  Soc  Japan,  51 : 450-57 
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the  transpose  of  the  noise  vector,  and  e.  is  a  normally  distributed  random 

variable  with  mean  zero  and  variance  o2,  then  the  regression  weights 
determined  from  the  data  are 

6  *  [xVxr  [xV’v]  ’  <3, 

where  -1  denotes  matrix  inversion  and  the  objective  analysis  may  be 
accomplished  by 


See  Draper  and  Smith6  for  details. 

The  kth,  1th  element  of  the  matrix  [ X tV ~ 1 X 1  is  computed  by 


i 


N  N 

l  I 

i  j 


X..X  . ,  V  7 
ik  jl  ij 


(5) 


and  the  kth  element  of  the  vector 


[xVW]  is  computed  by 


N  N 

I  l  X..Y.V7 
i  j  u  J  1J  ' 

(6) 

The  covariance  structure  function,  cov(XX),  is  determined  from  the  histor¬ 
ical  data  set  for  all  possible  kth,  1th  elements  of  [XtV'1X]  and  for  all 


6N.  R.  Draper  and  H.  Smith,  1966,  Applied  Regression  Analysis, 
Wiley  Press,  401  pp 
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possible  kth  elements  of  [XtV-1Y],  and  is  expressed  in  terms  of  a  mathe¬ 
matical  model.  The  signal  plus  noise  covariance  can  be  computed  from  this 
modeled  structure  for  any  kth,  1th  element  required,  within  the  limits  of 
the  model  itself. 


At  =  1 2  ~  1 1 

Xi  =  position  (E-W)  of  variable  1 
X2  =  position  (E-W)  of  variable  2 
Y:  =  position  (N-S)  of  variable  1 
Y2  =  position  (N-S)  of  variable  2 
ti  =  observation  time  of  variable  1 
t2  =  observation  time  of  variable  2 

The  decision  variables  (NLP  estimated  parameters)  are  as  follows: 

A  =  modeled  lag  zero  correlation  coefficient. 

a  ,  a  ,  c.  =  modeled  measures  of  storm  size  in  space  and  time, 
x  y  l 

a  =  a  measure  of  ellipticity. 

C  and  C  =  modeled  speeds  toward  the  east  and  north  (positive), 
x  y 

respectively  (figure  3). 

Crawford's7  NLP  algorithm  minimizes  an  objective  function,  subject  to 
mathematical  constraints.  The  algorithm  adjusts  the  values  of  the  vari¬ 
ables  until  any  minor  change  in  these  variables  results  in  an  improvement 
of  the  objective  function  of  less  than  some  small  preset  value.  The 
algorithm  thus  converges  upon  values  of  the  variables  which  are  as  good 
as  or  better  than  any  nearby  values. 

In  the  above  formulation,  variables  1  and  2  were  both  IR  or  visible,  thus 
producing  an  autocovariance  function. 


SYNOPTIC  SITUATION 

The  first  case  used  to  test  the  objective  analysis  procedures  occurred  on 
20-21  May  1977  over  the  south-central  United  States  (see  appendix  for  map 
series).  The  synoptic  situation  on  19  May  at  1200Z  showed  a  high  pressure 
system  centered  over  Kentucky  with  a  low  pressure  over  North  Dakota.  A 


7K.  C.  Crawford,  1977,  "The  Design  of  a  Multivariate  Mesoscale  Field 
Experiment,"  PhD  Dissertation,  University  of  Oklahoma,  Department  of 
Meteorology 
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Figure  3:  Cloud  top  characteristics  for  one  portion  of  the  storm,  lifetime 
deduced  from  the  modelled  covariance  function.  a=major  axis, 
b= minor  axis,  e  ^orientation,  c  x  =  eastward  speed  of  storm, 
c  y  =  southward  speed  of  storm,  *  =  direction  of  movement 
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cold  front  extended  from  this  low  pressure  system  southwestward  through 
South  Dakota,  Nebraska,  western  Kansas,  and  Oklahoma  into  the  Texas  Panhandl 
Flow  in  advance  of  the  front  was  southerly,  bringing  warm  moist  air  from  the 
Gulf  of  Mexico  in  the  Great  Plains.  Surface  temperatures  in  advance  of  the 
front  ranged  from  the  lower  70 ’ s  in  North  Texas  to  the  lower  60’ s  in 
North  Dakota  with  dew  point  temperatures  ranging  from  the  mid  60' s  to  the 
mid  50' s. 

Upper  airflow  at  700  mb  on  the  19th  was  characterized  by  a  ridge  over  the 
eastern  portion  of  the  United  States,  while  a  trough  spanned  the  area  from 
Idaho  through  eastern  Nevada  into  southern  California. 

Twenty-four  hours  later  the  surface  high  pressure  system  had  moved  over 
Pennsylvania,  and  the  cold  front  reached  from  northwest  Wisconsin  to  central 
New  Mexico.  Warm  moist  Gulf  air  continued  to  flow  into  regions  in  advance 
of  the  cold  front.  The  upper  level  trough  and  ridge  positions  had  moved 
slowly  eastward. 

On  the  mesoscale,  a  line  of  thunderstorms  had  developed  by  1700Z  in  the 
Texas  Panhandle  with  the  line  moving  southeastward  and  the  cells  in  the 
line  moving  northeastward .  Later  in  the  day  multiple  lines  had  formed 
and  thunderstorm  tops  to  14  km  were  reported  (see  appendix  for  radar 
echo  maps).  In  addition,  rainshowers  developed  around  1630Z  in  central 
Texas,  and  these  two  areas  me'  ged  by  1930Z  over  the  Red  River.  The  greatest 
reported  hourly  rainfall  amount  occurred  between  0000Z  and  0100Z  on  21  May 
at  Geary,  Oklahoma  (2.00  in.).  Hourly  rainfall  amounts  in  excess  of  1  in. 
were  recorded  from  many  areas  in  the  Texas  Panhandle  and  western  Oklahoma 
between  1800Z  on  20  May  and  0100Z  on  21  May.  Total  storm  rainfall  (see 
appendix  maps)  pattern  was  oriented  southwest-northeast  with  the  largest 
amount  (4.3  in.)  occurring  at  McLean,  Texas.  The  3-in.  isohyet  ran  from 
near  Kress,  Texas,  northeastward  to  Sun  City,  Kansas. 


AUTOCOVARIANCE  AND  CROSS-COVARIANCE  FUNCTION  RESULTS 

The  satellite  data  for  the  20  May  1977  storm  produced  the  values  shown 
in  part  A  of  tables  1  and  2.  Storm  features  (part  B)  derived  from  the 
values  in  part  A  were  computed  in  the  following  manner.  The  orientation 
(e)  is  given  by 

6  =  (1/2)  ARCTAN  [2aox  oy/(o2  -  o2)]  (8) 


The  major  axis  (a)  is  computed  from  equation  (9) 


a2  =  a2  cos2  e  +  a2  sin2  e  +  2a a  o  sin  e  cos  e  .  (9) 

x  y  x  y 
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TABLE  2.  VISIBLE  DATA  20  MAY  1977 


0.34  0.58  2100  -  2145  15  280  140  123  2.00  224  64  0.51 


The  minor  axis(b)  is  given  by 


b2  =  o2  sin2  0  +  o2  cos2  0  -  2ao  o 
x  y  x  y 

The  area  is  computed  from  equation  (11) 


sin  e  cos  0  . 


A  =  7T  (a )  ( b)  ; 


(10) 


(ID 


and  the  shape  by 


=  a/b  ; 


(12) 


and  the  signal  to  noise  ratio  (SNR)  can  be  inferred  from  the  modeled  lag 
zero  correlation  coefficient  as  follows  (equation  (13)) 


SNR 


1  -  pz 


(13) 


Also,  C  is  the  speed  of  movement  of  the  storm  and  <+>  is  the  direction  of 
movement. 


The  derived  storm  characteristics  are  conservative  and  stable  over  the 
1800Z  to  2145Z  time  period.  The  change  in  orientation  shown  at  2000Z  was 
due  to  the  merging  of  a  small  cell  with  the  convective  complex  (see  appen¬ 
dix  for  radar  coverage).  The  colder  IR  temperatures  became  numerous  in 
the  southern  portion  of  the  storm  during  the  merger,  thus  forcing  the 
orientation  angle  around  to  a  more  easterly  direction;  however,  the  direc¬ 
tion  of  movement  (<f>)  of  the  total  storm  remained  about  the  same.  Figures 
4  and  5  show  a  plot  of  some  of  the  significant  parameters  of  the  convective 
complex.  Weather  radar  reports  from  Lubbock,  Texas,  indicated  average 
cell  movements  of  30  to  35  knots  (55  to  64  km/hr)  with  a  direction  of  230 
at  1800Z.  Maximum  observed  cell  speed  by  weather  radar  was  50  knots 
(91  km/hr)  from  230.  Later  in  the  day  (2030Z),  the  Oklahoma  City  weather 
radar  reported  cells  moving  at  about  25  to  30  knots  (46  to  55  km/hr)  from 
230  to  240  degrees  azimuth.  A  maximum  speed  of  35  knots  with  an  azimuth 
of  230  was  also  recorded. 


The  computed  storm  orientation  (e)  agrees  well  with  the  surface  rainfall 
pattern  up  until  time  of  cell  merger.  .'(See  appendix  for  rainfall  analysis.) 
The  area  shown  in  tables  1  and  2  indicates  that  the  visible  cloud  top  is 
somewhat  smaller  and  oriented  a  little  differently  than  the  IR  cloud  top. 

The  overshooting  tops  seen  in  the  visible  data  (figure  6)  are  pronounced 
and  tend  to  line  up  along  the  direction  of  motion.  These  features  tend 
to  make  the  covariance  algorithm  deduce  a  smaller  system. 
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STORM  ORIENTATION  (0)  STORM  SPEED  (Km/Hr) 

Figure  4.  Significant  parameters  of  convective  complex. 
IR  threshold  data  =  -30°C. 


ORIENTATION(DEG) 


DIRECTION(DEG) 


Figure  5 


Significant  parameters  of  convective  complex, 
Visible  data  20  May  77. 


The  SNR  increased  and  then  decreased  with  time.  The  convective  complex 
evolved  during  the  period  of  analysis  from  a  smooth  and  well-defined 
ellipse  up  until  the  time  of  cell  merger.  Then  the  system  became  a  rather 
flat  and  ragged  composite  of  anvils  produced  by  several  embedded  cells 
(figure  7).  The  SNR  could  be  anticipated  to  drop  to  very  small  values 
as  the  system  ceased  to  produce  new  active  rain  cells. 

Several  inferences  can  be  made  from  the  results  of  this  study,  namely: 

(1)  cloud  top  morphology  of  convective  complexes  as  deduced  by  using  the 
space-time  covariance  analysis  technique  is  sufficiently  conservative  to 
suggest  that  predictions  made  for  periods  of  about  1  to  2  hours  could 
show  significant  skill;  (2)  these  predictions  can  be  objective  and  auto¬ 
mated;  ana  (3)  rainfall  location,  movement,  and  intensity  can  be  associated 
with  cloud  top  features  of  such  weather  systems. 

All  the  statistics  required  to  perform  the  analysis  of  the  GOES  data  are 
contained  in  the  cross-covariance  functions  derived  from  the  data  sets. 

The  first  step  in  modeling  these  functions  is  to  obtain  "raw,"  discrete 
covariance  estimates  by  lagging  each  field  with  respect  to  itself  in  space 
and  time.  To  cover  the  entire  storm  and  yet  preserve  an  adequate  space- 
time  resolution,  the  raw  covariance  matrices  were  of  order  9  by  9  in  space 
and  4  in  time.  The  category  sizes  or  resolution  were  70  km  east-west, 

70  km  north-south  and  15  minutes  in  time.  The  product  pairs  and  the  raw 
covariance  matrices  for  selected  times  are  presented  in  the  appendix. 

The  raw  covariance  matrices  are  normalized  to  give  raw  correlation  matrices. 
Equation  (7)  was  fitted  to  the  raw  correlation  matrices  by  using  a  non¬ 
linear  program  (NLP)  algorithm  developed  by  Crawford.7 


Cov  (XjX2)  =  A  *  exp  j  [-1/2(1  -  a2)]  [(&X/ox)2 

-  2a(AX/a  )  (AY/a  )  +  (aY/o  )21 

*  y  y  J 

-  (1/2)  (At/at)2J*  COS  |  (it/2)  [(AX/ax)2 

-  2a(AX/ax)  (aY/o^)  +  (AY/o^)2  +  (At/c^)2]1/2  J  (7) 


7K.  C.  Crawford,  1977,  "The  Design  of  a  Multivariate  Mesoscale  Field 
Experiment,"  PhD  Dissertation,  University  of  Oklahoma,  Department  of 
Meteorology 


where 


AX  =  (x2  -  Xj)  -  cxt  -  xo 
AY  =  (Y2  -  YO  -  Cyt  -  Yo 


RAIN  GAGE  -  SATELLITE  DATA  COMPARISON 

The  procedures  in  the  previous  action  showed  that  an  ellipse  can  be  fitted  to 
a  convective  complex  and  useful  parameters  can  be  derived  from  the  auto¬ 
covariance  and  cross-covariance  functions.  However,  if  an  analysis  system 
is  to  operate  in  a  real-time  mode,  the  above  procedures  require  too  much 
computer  time  for  operational  purposes. 

An  ellipse  fitting  routine*  was  developed  for  the  AMAS  image  processor  which 
allows  an  ellipse,  similar  in  configuration  to  the  ones  in  the  previous 
section,  to  be  fitted  to  the  convective  complex  temperature  patterns  in  a 
few  seconds  of  computer  time.  Successive  images  can  be  averaged  or  dis¬ 
played  individually  as  the  analyst  desires.  This  system  also  allows  image 
overlay  of  other  data.  Using  the  navigation  system,  one  can  take  each 
line-element  location  on  the  image  and  transform  it  into  latitude  and 
longitude;  one  can  transform  earth-located  information  into  the  satellite 
coordinate  system  or  vice-versa  for  proper  comparison  of  the  data.  This 
system  permits  new  insight  into  the  physical  mechanisms  occurring  during 
the  rainfall  process  and  demonstrates  where  efforts  should  be  made  in 
further  studies. 

Rainfall  data  for  the  20  May  1977  period  was  available  at  hourly  intervals 
from  some  stations  and  on  a  daily  total  basis  from  other  stations.  Figures 
8,  9,  10,  and  11  compare  the  hourly  rainfall  totals  for  1900Z,  2000Z,  2100Z, 
and  2200Z  with  the  fitted  ellipses  which  correspond  to  the  -63°C,  -65°C, 
and  -67°C  IR  temperatures  averaged  for  the  half-hour  and  hourly  values  end¬ 
ing  on  the  hour.  Since  the  available  satellite  data  did  not  span  a  24-hour 
period,  no  comparison  was  made  with  the  total  rainfall  pattern.  The  greatest 
amounts  of  recorded  rainfall  fell  within  the  -67°C  fitted  ellipse.  All 
but  two  rainfall  events  greater  than  0.50  in.  fell  within  the  -67°C  ellipse. 
However  one  station  (within  the  -67°C  ellipse)  did  not  report  any  rain  dur¬ 
ing  the  4-hr  time  period.  The  -63°C  fitted  ellipse  encompassed  most  of  the 
significant  (>  0.50  in. /hr)  rainfall  events.  Only  one  event  fell  outside 
this  ellipse.  More  rain  usually  fell  after  mergers  of  -67°C  contour  areas. 


*J.  Marvin,  Physical  Science  Laboratory,  New  Mexico  State  University, 
personal  communication 
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Figure  8.  Hourly  rainfall  rate,  1900Z,  20  May  77,  compared  with 
fitted  ellipses. 


Figures  12  and  13  show  the  shape,  position,  and  pixel  areas  of  the  -67°C 
contours  for  15-min  intervals  from  1800Z  to  2145Z.  If  the  centroid  of  the 
-67°C  is  plotted  with  respect  to  time,  speed  of  movement  (43  to  54  km/hr) 
of  the  centroid  is  very  close  to  the  radar  detected  speed  of  movement  of 
the  thunderstorm  cells. 

From  this  study,  it  appears  that  certain  ellipses  fitted  to  IR  temperature 
levels  can  be  assigned  rainfall  amounts.  Also,  more  rain  can  be  assigned 
to  the  downwind  half  of  the  ellipse  than  to  the  upwind  half.  However,  if 
this  scheme  is  to  be  automated,  the  subjective  judaments  of  when  area  mer¬ 
gers  occur  and  location  of  downwind  areas  of  ellipses  must  be  eliminated. 
Therefore,  the  following  table  is  suggested  for  use  in  solving  for  rainfall 
values  to  be  used  in  trafficabil ity  calculations. 


TABLE  3.  RAINFALL  ESTIMATIONS  FOR  VARIOUS  ELLIPSES 


Ellipse  (°C) 

Rainfall  Range  (in. /hr) 

-67 

0.50  to  1.50 

-65 

0.25  to  0.50 

-63 

Trace  to  0.25 

Although  these  values  probably  overestimate  the  rainfall  in  some  areas  and 
underestimate  it  in  others,  they  do  yield  reasonable  values  without  sub¬ 
jective  judgments.  A  problem  area  in  this  type  of  estimation  is  rainfall 
efficiency.  The  base  of  the  clouds  is  often  several  hundred  meters  above 
the  surface  and  subcloud  evaporation  of  rainfall  is  often  quite  high,  which 
means  that  these  same  cloud  top  temperatures  could  yield  different  rainfall 
values  in  other  areas  and  times  of  the  year.  Possibly  cl imatoloqies  of 
rainfall  estimates  will  have  to  be  made  according  to  time  of  year,  time  of 
day,  and  location  before  accuracy  can  be  improved.  Additional  case  studies 
must  be  done  to  accomplish  the  above  task. 

CONCLUSIONS 

The  IR  data  from  the  cloud  tops  of  convective  complexes  relate  well  to  rain¬ 
fall  at  the  surface.  Maximum  rainfall  occurs  when  cloud  top  temperatures 
are  -67°C  or  colder. 

The  centroid  of  the  fitted  ellipse  of  the  -67°C  contour  area,  when  plotted 
against  time,  appears  to  agree  well  with  the  speed  of  movement  of  cells 
tracked  by  weather  radar. 

Cloud  top  morphology  can  be  detected  by  using  the  AMAS,  and  ellipses  (fitted 
to  the  data)  can  be  drawn  directly  on  the  display  screen. 


CO'  KN 


Finure  12.  Shape,  position,  and  pixel  areas  of  the  -67°C  contours  for  15-min  intervals. 

Numbers  under  the  contour  areas  indicate  the  number  of  satellite  pixels  withii 
the  contour.  Mergers  of  areas  are  indicated  by  letter  desipnator  chance  from 
one  time  period  to  another. 


Convective  complex  movements  and  areal  expansion  rates  can  be  predicted  by 
using  past  positions  of  the  centroid  and  areas  of  a  fitted  ellipse  and 
extrapolation  techniques. 
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Figure  A-3.  Surface  at  2000Z,  20  May  77. 
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Fioure  A-4.  Surface  at  2300Z,  20  May  77. 
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